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ABSTRACT 

We address whether possible scale-dependent deviations from gaussianity in the pri- 
mordial density field that are consistent with the cosmic microwave background obser- 
vations could explain the apparent excess of early cluster formation at high rcdshift. 
Using two phenomenological non-gaussian models we find that at fixed normalisation 
to the observed local abundance of massive clusters, the protoclusters observed at z^4 
are significantly more likely to develop in strongly non-gaussian models than in the 
gaussian paradigm. We compute the relative z < 1 evolution of X-ray cluster counts 
in the non-gaussian case with respect to the gaussian expectation, and the relative 
excess contribution to the CMB power spectrum due to the integrated thermal SZ 
effect. We find that both the observed hints of an unexpectedly slow evolution in the 
X-ray counts and the excess power at high £ that may have been observed by CMB 
interferometers can also be reproduced in our non-gaussian simulations. 
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1 INTRODUCTION 

Non-gaussianities in the primordial perturbations of the cos- 
mic density field are inherently difficult to generate in simple 
versions of inflation, and most inflationary models predict 
negligible amounts of relic non-gaussianity. Multiple field 
inflation or more speculative models such as those assuming 
non-homogeneous reheating processes or those seeded with 
topological defects can however readily rise to the challenge 
of spawning deviations from a gaussian initial density field. 
Of course, the amount of allowed primordial non-gaussianity 
on CMB scales has come u nder increasing pressu re from the 
recent results of WMAP llKomatsu et alJl2003ft . and con- 
straints are expected to tighten with upcoming data releases. 

However current constraints from all-sky CMB maps 
only transfer to constraints at cluster or galaxy scales in the 
case of scale-invariant non-gaussianity. If non-gaussianity is 
strongly scale-dependent and arises only at cluster scales, it 
can simultaneously fit the WMAP bounds yet affect large- 
scale structure so as to induce predictable deviations from 
the gaussian paradigm. The caveat here is the "fine-tuning" 
of the physical scale at which non-gaussianity appears, which 
has to fall between I ~ 300 and I ~ 2500. (Models including 
defects such as linear cosmic strings as additional seeds for 
structure formation would naturally have such a scale, see 
lAvelino et alj Fl998) . Our approach here is phenomenolog- 
ical, and we consider two hypothetical models which have 
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such scale-dependent non-gaussianity, and, incidentally, are 
also to some extent motivated by theoretical considerations. 

From the point of view of structure formation at low 
redshift and relatively small scales, issues for the gaussian 
ACDM paradigm include the early formation of protoclus- 
ters up to z ~ 4 as traced by radio galaxies surrounded by 
Ly-a emitters and Lyman-break galaxies, the observed ap- 
parently slow evolution of the comoving density of bright 
X-ray clusters, and the possible excess of angular power in 
CMB temperature fluctuations seen by CBI and ACBAR 
on cluster scales that most likely is due to the thermal SZ 
effect. It is well known that non-gaussianity may be able to 
naturally explain possible features associated with an excess 
(relative to Gaussian predictions) in the frequency at high 
redshift of galaxy protoclusters and clusters, by reasoning 
as follows. Rare peaks in the primordial density field de- 
fine the locations where massive clusters will form, and the 
higher the peak the earlier the collapse to a virialized sys- 
tem. If the probability density function (PDF) of the den- 
sity field shows an excess of high peaks compared to the 
gaussian case, (proto-) clusters of given mass form earlier, 
and systems above a given mass will be more abundant at 
fixed redshift. Note t hat more speculat i ve top ics such as the 
periodicity noted bv iBroadhurst et al.l ll990|), the p r esence 
of large regions underdense in galaxies iFrith et al.l (j'2003) 
or the high correlation length of very massive X-ray clus- 
ters could also hint at primordial non-gaussianity. Note also 
that scale-dependent no n-ga ussianity has been advocated by 
lAvelino fc Liddld j2004) and lChen etafl j2003ft as a possible 
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way to reconcile the high optical depth to reionization with 
realistic efficiencies for the formation of the first stars, and 
that its possible impact on th e inner profiles of dark m atter 
haloes has been addressed bv lAvila-Reese et al.l J2003T) . 

Here, we use simulations of large-scale structure to 
check the impact of strongly positively and negatively 
skewed non-gaussianity on the above observables. We evolve 
two primordial cosmic density fields with a scale-dependent 
non-gaussianity which is at least of the same order as the 
gaussian component at the scales probed by the simulation 
box. We employ the concordance cosmological parameters 
and normalise the models so that they give a similar abun- 
dance of massive clusters at z = 0. The two non-gaussian 
models differ by the sign of their primord ial skewness: we 
choose a model with primordial voids dGrifnths et aljl2003t) 
to realize negative skewness, and a model with a chi-square 
PDF of the density field to have positive skewness ( this case 
been worked out in the scale-independent case bv iPeeblesI 
1999). We compare to the gaussian ACDM case they bracket 
from the point of view of cluster formation. 

The paper is organised as follows. In Section |5] we re- 
view recent observational hints of non-gaussianity on cluster 
scales. Section [3] describes the setup of the simulations and 
compares to observations. We conclude in Section ^] 



2 INDICATIVE OBSERVATIONS 

We describe here the three observational hints which we can 
probe using simulations without recourse to any model for 
galaxy formation. 

2.1 High redshift assembling structures 

IWillickl i200(J) was among the first to set constraints on pos- 
sible primordial non-gaussianity using the presence of high- 
redshift massive systems. Even if some of the mass estimates 
he used fo r the z = 0.83 cluster MS 1054-03 were proba - 
bly biased jHoekstra et alJl200d:lNeumann fc Arnaudl200oT) 
because the system might not be dynamically relaxed, the 
precedent was set. 

On t he theoretica l side , [ Robinson et af| J2000I) 
(see also iKovama et alj 1 19991 : iRobinson fc Bakerl l2000t 
lMatarre"seet"^TT200oF have shown that given fi ma tter,o, 
the power spectrum and its normalisation to the cluster 
abundance today, the mere presence of massive high redshift 
clusters provided useful constraints on relative amounts 
of generic non-gaussian contributions to the primordial 
density field. 

More recently, iMilev et aD <2004l) (see also iKurk et all 

2003) report the detection of 21 Ly-a emitters assembling 
around a radio galaxy at 2 ~ 4.1. The full structure has 
a proper spatial extent of 3 t o 5 M pc. Even if the system 
is not virialized, IMilev et alJ i2004h could compute its for- 
mal velocity dispersion from the Ly-a systems and found it 
to be of order 325 km s _1 . The enclosed mass is compara- 
ble to that of local galaxy clusters. At z ~ 2.1. IKurk et alJ 
( 2003) find a close pair of two such Ly-a systems; when they 
are considered seperately, they have velocity dispersions of, 
respectively, ~ 200 and ~ 400 km s _1 . When they are con- 
sidered a single group, the estimated velocity dispersion is 
1000 km s" 1 . 



2.2 Evolution of the X-ray cluster luminosity 
function 

The X-ray window is ideal for studying the evolution of mas- 
sive clusters because the selection function effects can be 
much better understood than in the optical. No direct mass 
measurement is needed and one can use the observed rela- 
tions between the X-ray luminosity and temperature (L — T) 
and between the temperature and mass (M — T) to select 
massive members from a catalogue of X-ray clusters. In- 
evitably however, results on the evolution of the abundance 
of X-ray clusters will depend on the adopted scaling rela- 
tions. 

Evolution of the abundance of X-ray luminous clus- 
ters has become one of the major tools for constraining the 
present d ay cosmological m atte r density S7 m atter.o = 1 (for 
instance. ISadat et alJll998l see iBoreani etalll200lj for ev- 
idence for a low-density universe). However, at fixed cos- 
mological parameters, fixed scaling relations and assuming 
normalisation to the local cluster abundance, increasing the 
amount of skewness in the primordial density field changes 
the z <J 1 evolution in the abundance of massive clusters. 

Deep ROSAT, Chandra and XMM X-ray observations 
have revealed a significan t number of distant luminous clus- 
ters llRosati et alj |2000). At redshifts z <J 0.8, statistics 
of X-ray clusters have become sufficiently good to pro- 
vide realistic e s timat es of their cosmological a bund ance . 
IVikhlinin etaD (1200311 : IVoevodkin fc Vikhlininl <2004h de- 
rive cosmological constraints from the evolution of the clus- 
ter baryon mass fraction which they relate to the total clus- 
ter mass, thereby skipping the usual uncertainties in the 
scaling relations. Using the 160 deg 2 ROSAT survey, they 
obtain results in line with the evolution of the mass function 
expected in a concordance gaussian CDM cosmology as they 
find moderately negative evolution of the cluster abundance 
to z ~ 0.6. 

Recently. IVauclair et alJ [20031 prefer an f2 mat tcr,o = 1 
universe from the reduction of the same data, but caution 
that a variation in the scaling of the M — T relation with red- 
shift might bias the analysis. In this work we will rederive 
the evolution of the density of lum inous clusters fr o m the 
~200-member cluster catalogue of Vikhlinin^^lJ_J^98jj 
(updated using spectroscopic redshifts in lMullis et alJl2003i : 
we will refer to this catalog as VM) and we will compare it 
with predictions from gaussian and non-gaussian models, at 
fixed scaling relations. 

2.3 Small-scale CMB power spectrum 

iMaiumdar fc Subrahmanvanl i2000l) used the upper bounds 
on the rms CMB temperature fluctuations observed at ar- 
cmin scales by ATCA as constraints for structure formation 
models. They ruled out the possibility of a high f2 ma ttcr,o 
assuming COBE normalisation and gaussian fluctuations 
of the primordial density field. More recently, high reso- 
lution CMB experiments suc h as CBI iMason et al] I2OO3T) 
or ACBAR llKuo et alJl2004l) find an excess of small-scale 
(I ~ 3000 or 8 « few arcmin) power with respect to the 
expectations from primary CMB fluctuations. 

Non-gaussian models deviate more from gaussian mod- 
els at high redshift (at fixed present day normalisation). 
Because the SZ temperature decrement of the CMB is in- 
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dependent of redshift, the angular power spectrum of the 
thermal SZ effect is sensitive to the high redshift abundance 
of clusters. This turns the integrated SZ power spectrum into 
another efficient means of distinguishing gaussian from non- 
gaussian models, in addition to the evolution of the abun- 
dance of massive clusters. 

Given the known local abundance of clusters, gaussian 
models require low values of the normalisation (as ~ 0.8 if 
fi mat ter,o = 0.3) when they are fitted to the data. It is not 
clear yet if the observed excess is due to systematic effects or 
if it is a real signature coming from compact sources (point 
sources and/or galaxy clusters). In case this excess in the 
CMB power spectrum at I > 3000 is confirmed by future 
experiments to be due to galaxy clusters, it would imply 
a high normalisation (as > 1 for fi ma ttcr,o ~ 0.3) in the 
gaussian case. Alternatively, it would point towards a com- 
bination of the M — T and r coro — M relations (detailed in 
paragraph s. 3. 3lbelow) which would be rather extreme (see 
for example lMaiumdarll200ll) . A more natural explanation 
however is an over-abundant population of clusters at high 
redshift, in turn a possible consequence of non-gaussianity 
of the primordial density field. We will consider this possi- 
bility in line with the evolution of the scaling relations that 
we assume to compute the X-ray cluster evolution. 



3 THE CASE FOR NON-GAUSSIANITY 

We first describe the simulations and characterize the mod- 
els with simple statistics. Then, we directly measure in the 
simulations the three observables of the previous Section. 

3.1 Simulations setup 

We consider two non-gaussian models: the fi rst, called V is 
a model with compensat ed primordial voids iGriffiths et alJ 
120031: iMathis et al]l2003h . the second (hereafter model C 2 ) 
is a model where the PDF of the primordial field density 
is the square of a gaussian. In theory the models could be 
associated with, respectively, a possible first-order transi- 
tion occurring during inflation, and to versions of multiple 
scalar field inflation with isocurvature-type initial perturba- 
tions in a massive scalar field playing the role of CDM (suc h 
an ICDM model was worked out in detail bv lPeebleslfl99g() . 
While the first-order phase transition model yielding primor- 
dial voids would be truly gaussian on CMB scales, this is not 
the case for the ICDM model. Tuning the ICDM model ini- 
tially suggested so that it has the correct scale-dependent 
non-gaussianity is beyond the scope of this paper. However 
in conformity with our phenomenological approach, we will 
refrain in the following from linking these V and C 2 reali- 
sations of the initial density field to a specific inflationary 
scenario. Rather, we simply take them as possible represen- 
tations of the primordial density field on cluster scales that 
could arise in generic inflationary models. For reference, we 
also consider a gaussian model G. 

All three models share the same ACDM cosmological 
parameters S7cdm,o = 0.3, Oa.o = 0.7 and h = 0.7. We im- 
pose the same initial matter power spectrum in the V and 
G models: n s = l with an adiabatic CDM transfer function. 
In the C 2 mode l, the initial power spectrum is that used 
bv lPeebleslll999l with n s =-l.8 with an isocurvature CDM 
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Figure 1. Mass function of dark matter haloes at z = for our three 
phenomenologi cal models. Diamond s with error bars (resp . triangles) 
are the data of lBahcall et alJ i2003ft jVikhlinin et all2003ft . 
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Figure 2. Mass function of dark matter haloes at z = 4 for our three 
phenomenological models. 



transfer function. We ensure that the z = cluster mass 
functions of the non-gaussian models match reasonably well 
those of the gaussian ACDM model with as = 0.9. This is 
shown in Fig. Q where we compare the mass functions in 
our three si mulated model s to th e local optical cluster mass 
function of iBahca ll et al.| i2003j) and to the X-ray cluster 
mass function of lVikhlinin et aTT(l20 03') . Differences with re- 
spect to the gaussian case are limited to 25 percent up to 
1.6 x 10 14 /i _1 M e cluster masses. At masses ~ 10 15 /i _1 M Q , 
differences reflect cosmic scatter. In practice, the match to 
the ACDM G model was achieved in the C 2 model by nor- 
malising the initial power computed from the realized field 
so that linear extrapolation yields as = 0.65 today. In V, 
the underlying gaussian field before "addition" of the voids 
has a normalisation similar to that of G. Note that the pa- 
rameters of the distribution of void sizes in the V model 
was set so as to matc h at present day the large voids seen in 
local galaxy surveys fel-Ad fc Piranll2000tlHovle fc Voeelevl 
2003). As such, V would also be a natural solution to the 
presence of large regions in the near by Universe w hich are 
completely empty of bright galaxies dPeeblesll200"iTl . 
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The simulations start at high redshift z- m it = 300 and 
follow collisionless structure formation using 128 3 particles 
in a 100 h^ 1 Mpc box with periodic boundary conditions. 
We employ the publically available code N-body code GAD- 
GET without hydrodynamics 1 and set the equivalent Plum- 
mer softening length to one tenth of the mean interparticle 
separation. The mass of a dark matter particle is 3.94 x 10 10 
h' 1 Mp>. 



3.2 Non-gaussian signatures in the simulations 

This specific choice of the strongly non-gaussian models V 
and C 2 is attractive in the sense that the skewness of the pri- 
mordial density field approximately brackets that of a gaus- 
sian field. This is illustrated in Fig.|3|where we plot the PDF 
of the rescaled overdensity field S/a directly measured it in 
the three simulations at z — 50 over 8 ft -1 Mpc scales. Dot- 
ted, solid and dashed lines respectively correspond to the V, 
G and C 2 cases. For comparison, a gaussian has been over- 
plotted in dash-dotted line. While the PDF of the density 
of the V simulation is strongly negatively skewed due to the 
presence of the voids, that of the C 2 simulation is positively 
skewed as expected. (Note that the C 2 and G simulations 
are still within the "linear regime" over 8 h~ x Mpc scales 
by z = 50, while the expanding primordial voids in the V 
simulation have already fully emptied regions as large as 8 
h,- 1 Mpc.) 

Such skewed primordial PDFs directly impact the z — 
4 halo mass function as shown in Fig. [5] There, we find 
early structure formation to be a common feature of the 
C 2 and V models. In the V model, this is the result of two 
effects: fragmentation of the shell surrounding the voids, and 
non-linear motions in the background which are induced by 
the void underdensities. In the C 2 model, it is the effect of 
the relative primordial overabundance of > 3u peaks with 
respect to the gaussian case. 

The effect of the three different initial conditions is also 
clearly mirrored, although differently, in the underdensity 
proba bility function (UPF) UPF <(5 (r) JWeinberg fc Cold 
1992). We define the UPF as the probability that the matter 
underdensity 8 in a randomly placed spher e of comoving ra- 
dius r has a value less that -0.8, following iHovle fc Vogelevl 
( 2003). Fig.^Jshows the underdensity probability function at 
z = for the gaussian G, the C 2 and the V models (solid, 
dashed and dotted lines respectively). Obviously, the void 
model has a larger probability of underdensities compared 
to the gaussian case, while the C 2 model has very few re- 
gions with a radius larger than 1 h~ Mpc where <5 5j —0.8, a 
consequence of the \ 2 distribution of the PDF of the primor- 
dial field. A comparison to results from local galaxy surveys 
is nevertheless out of the present scope: this would require 
a model for galaxy bias, which plays a critical role in void 
statistics. 



3.3 Comparison to observations 

We first compute statistics relevant to the protoclusters ob- 
served at z ~ 4. Then, we choose two models for the cluster 
scaling relations and employ them to compare the observed 
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Figure 3. PDF of the rescaled overdensity field S/a measured at 
z = 50 in the simulations, for the two non-gaussian and for the gaus- 
sian model. A gaussian distribution (dash-dotted line) has been over- 
plotted. Note the opposite skewness of the C 2 and V models. 




Figure 4. Underdensity probability function UPF<^(r) for S = —0.8 
measured at z = in our three models. Note the differences in the 
non-gaussian models V and C 2 (dotted and dashed lines) compared to 
the gaussian case (solid line). The void model naturally has a higher 
UPF(r), while the C 2 model has very few underdense regions larger 
than 1 h~ 1 Mpc in radius. 



evolution of the X-ray luminosity function. Finally, we es- 
timate the expected thermal SZ angular power on cluster 
scales. 



3. 3. 1 Finding high-redshift protoclusters 

Comp aring the "protocluster" observations of iMilev et alJ 
(2004) to simulations of structure formation requires suit- 
able statistics. We proceed as follows. At z = 4, we ran- 
domly throw a large number (30000) of spheres of proper 
radius r^i-ow in the simulation and compute for each the 
overdensity with respect to the background and the one- 
dimensional velocity dispersion a v of all the dark matter par- 
ticles with respect to the centre of mass. We then construct 
the PDF of a v gathering all the spheres o f adequate overden- 
sity. (This procedure was suggested by iRob mson fc Bakerl 
2000 in the case of the sampling of a density field.) Note 
that this method does not require the overdensities within 
the spheres to have virialized. 

Fig. Qj] shows the resulting PDF at z = 4. The upper 
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left panel corresponds to throwing spheres with r t hrow = 5 
h~ x Mpc comoving and gathers the PDF of a v measured 
in the spheres irrespective of their environment (overden- 
sity with respect to the background) . Th e 3 to 5 Mpc 
proper extent reported by iMilev et al] 12004 maps on to a 
5 to 9 h -1 Mpc comoving sphere radius. The upper right 
panel of the same figure corresponds to throwing spheres 
with r t hrow — 9 h -1 Mpc again selected in all environments. 
The lower row of Fig.|S]is an addition to the upper row, for 
spheres whose dark matter overdensity 5dm £ [1,3]. Dot- 
ted, solid and dashed li nes are for the V, G and C 2 simu- 
lations respectively. The lMilev et alJ (120041) 325 km s" 1 ve- 
locity dispersion measure for the system of 21 Ly-a emitters 
at z = 4.1 is shown as the vertical dash-dotted line. Fig. El 
repeats the same statistics at z ~ 2, where other such pro- 
toclusters have been observed . The 200 and 400 km s _1 ve- 
locity dispersion estimates of iKurk et alJ (|2003f) are shown 
with vertical dash-dotted lines. (For the purpose of compar- 
ison we have kept the 5 and 9 ft -1 Mpc comoving radii at 
z = 2.) 

It is clear from these figures that at z = 4, there is no 5 
h~ x Mpc comoving radius sphere in the gaussian simulation 
G whose inner velocity dispersion is within 100 km s _1 of 
what is observed, even regardless of environment. The situa- 
tion is much more favourable for the other two non-gaussian 
models which we explore, where a non-zero PDF overlaps 
with the observations. The hypothesis here is clearly the 
absence of velocity bias between the simulated dark matter 
and the Ly-a emitters. For rthrow = 9 hT 1 Mpc, none of our 
simulations matches the observations. 

3.3.2 Simulating the abundance of luminous X-ray 
clusters 

We first select the most luminous clusters (Lx > L„) from 
the VM catalog and we build their density n(> L„) as a 
function of redshift: 

= (i) 

i 

where L* is the threshold luminosity and the sum is over 
all the clusters above L* in the redshift bin i (we take 
Az = 0.1). Vi is the search volume defined by the product of 
the volume in the shell of Az — 0.1 at redshift z times the ef- 
fective sky coverage f(L, z). For the 160 deg 2 ROSAT survey 
used here, f(L, z) was computed in lVikhlinin et aT] lll99cl) . 
In figure[7]we show results for a threshold L» = 2 x 10 43 /i -2 
erg s _1 cm -2 . Note that because we employ Ho = 100ft 
km s _1 Mpc -1 throughout, the original luminosities given 
for h — 0.5 in the VM catalog have to be rescaled. For ex- 
ample, 2 x 10 43 h~ 2 erg s 1 cm 2 correspond to an original 
0.8 x 10 44 erg s^ 1 cm" 2 for h = 0.5). 

Fig. [7| compares the observed n(> L t ) with our non- 
gaussian simulations (dotted and dashed lines are for V and 
C 2 respectively) and the gaussian G simulation (solid lines) . 
To make proper comparison between the models, the non- 
gaussian abundances have been slightly rescaled to have ex- 
actly the same z — abundance of L* clusters as in the 
G simulation. To take into account the uncertainties in the 
cluster scaling relations we have compared two different as- 
sumptions on the evolution of the scaling relations. Lower 
and upper lines are respectively for a no evolution model 




z 



Figure 7. Evolution of the density of luminous clusters wi th red shift. 
The sample of clust ers was taken from the lVikhlinin et alJ ll99Sl) and 
iMullis et alJ i2003t) catalog. Dotted and dashed lines show the ex- 
pected evolution for non-gaussian V and C 2 simulations respectively, 
while solid lines correspond to the gaussian case. In both cases, the 
upper line corresponds to a model where the L — M relation evolves as 
L ex M 2 (l + z) 3 while the lower line is for a model with no evolution 
(L oc M 2 ). 



L oc M 2 and for a model where L oc M 2 (l + z) 3 . The gaus- 
sian model could be compatible with the observed density 
of luminous clusters which is flat at least up to z — 0.6 only 
if there is a strong evolution in the L — M relation up to 
this redshi ft. It is interesting to compare our Fig . [7] with 
figure 9 in lMullis etaH (120041) . IMullis et alJ J2004I) present 
a curve similar to ours also using the 160 deg 2 sample, but 
they choose a higher L* = 10 45 /i^ 2 erg s _1 . They simi- 
larly find a decrease in the density of luminous clusters at 
z > 0.6. When they fit the observed density with an evolv- 
ing Schechter luminosity function normalised to match the 
z = 0.2 observed luminosity function of X-ray clusters, they 
recover this trend at z > 0.6. Their result is independent of 
any assumption (such as gaussianity) on the evolution of the 
mass function. In the non-gaussian case the L — M would 
need not to evolve or even to evolve negatively to explain 
the z — 0.6 drop. Whether this high-redshift drop is real or 
is due to the lack of statistics is an issue which has to be 
addressed with deeper and wider surveys. 

In fact, the L — M relation has been poorly determined 
so far, but we can use the recent attempts made to constrain 
the evolution of the L — T and M — T relation to constrain 
the evolution of the L — M relation. IVikhlinin et al.l <|2002T) 
constrain the evolution of the L — T re lation as L oc T°(l + 
z) 7 with a = 2.64 and 7 = 1.5 ± 0.3. iHolden et all (120021) 
find a softer value of 7 = 0.3 ±0.2. The BMW cluster survey 
jMoretti et all200ll) found evidence for a non-evolving L—T 
relation (or even a negative evolution). 
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Figure 5.2 = 4 PDF of the inner one-dimensional velocity dispersion measured from a large number of spheres thrown in the simulations. 
r throw = 5 h~ 1 Mpc is on the left, T t ^ TOVJ = 9 h~ 1 Mpc is on the right. Upper row: spheres in all environments, lower row: spheres where 
S G [1,3]. D otted, solid and da shed lines are for the V, G and C 2 simulations respectively. The vertical dash-dotted line is the system 
measured bv lMilev eta!] i2004h at z = 4.1. 



For the M — T relation. fvikhlinin et all §003) find T oc 
M x (l + z) i ' with V = 0.5±0.4. The spherical collapse model 
predicts tp = 1 for fimatter.o = 1. Combining the previous 
constraints one can derive that L oc M ax (l + z) 2-3 . 

iRosati et al] i2002l) gave a comprehensive review of the 
situation and suggest that a weak evolution in the L — T 
relation (7 < 1) is a reasonable assumption. If both the 
L — T and the M — T relation are non-evolving then the 
L — M relation must also be a non-evolving relation. With 
the current uncertainties in the observed cluster scaling re- 
lations we stress that it is unclear if the observed flatness 
in n(> L t ) is due to an evolving L — M or if it is due 
to an overdensity of high redshift clusters (when compared 
with gaussian models). Again, our purpose here is mainly to 
bring out the differences between gaussian and non-gaussian 
models at fixed scaling relations. 

We can use the following simple arguments to infer the 
redshift dependence of the L — M relation, if we assume the 
validity of a theoretical mass function model and place our- 
selves in the gaussian case. We use the fact that the luminos- 
ity function is compatible wit h no evolution up to redshift 
2 ~ 1 (for instance figure 9 in lRosati et, aljl2002f) . 

We simulated the luminosity function at different red- 
shifts for a flat cosmological model (as ~ 0.9, f2 ma ttor,o = 
0.3). The luminosity function can be derived from the mass 
function through the relation: 

dn dn dM , . 

d~L = lM~dL' ( ' 



Given the mass function (for example from the Press- 
Schechter formalism or its extensions to the non-gaussian 
cases) , the luminosity function only depends in the assumed 
L — M relation. Requiring that the luminosity function must 
be similar at low, intermediate and high redshifts we can 
constrain the redshift dependence of the L — M relation. In 
a simple test we have assumed the form: 



L = L M p (l + z)\ (3) 

changing j3 and £ and looking for the combination which 
yields a similar luminosity function at different redshifts (£ 
does not show a particular dependence with Lo). We have 
found that only models close to the relation £ = 1.5/3 — 1.1 
produce a non-evolving luminosity function. If j3 is derived 
from the observed T oc M 55 " 66 and L <x T 2 ' 5-3 then f3 w 
1.4 — 2. The corresponding £ for a non-evolving luminosity 
function reads £ w 1 — 2. This suggests that models with 
strong evolution (£ > 2) are not favoured by the observed 
non-evolving luminosity function. On the other hand, there 
must be some mild evolution in the L — M relations (£ J> 1) 
in order to reproduce the high redshift luminosity function. 
This reasoning holds for a gaussian primordial density field. 
In any case, if one finds using direct measurements of the 
mass and luminosity that £ < 1, it would then be difficult 
to explain the non-evolving luminosity function within the 
gaussian paradigm on cluster scales. 
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Figur e 6. Same a s Fig. [5| but at z=2. Vertical dash-dotted lines correspond to the ~ 200 and 400 km s velocity dispersion estimates of 
iKurk et al]<2003l) . 



3.3.3 Estimating the unresolved thermal SZ contribution 

Figure[S]illustrates the sensitivity of the SZ power spectrum 
at small angular scales to the abundance of the high red- 
shift population of clusters. In practice, the computation of 
the cluster SZ thermal power spectrum involves an integral 
of the mass function where each mass must be assigned a 
temperature. To compute the components Ct of the angu- 
lar po wer spectrum, we follow the formalism of iDieeo et alJ 
(2003[). We write the M — T scaling relation as: 

TocM x (l + z)*. (4) 

In this paragraph, the gas distribution inside the cluster 
needs modelling as well. For simplicity, we assume a /3-model 
with 13 = 2/3 and core radius r coro . We normalise the total 
emission to the total mass and the temperature. The angular 
power spectrum scales quadratically with the constant in 
front of equation 0| and we set it so that T — 8.5 keV for a 
M = 1O 15 /i _1 M cluster at z = 0. The results depend only 
mildly on \ and we fix it to x = 0.6. Our free parameters 
are therefore the redshift dependence of and the mass- 
and redshift-dependence of the core radius r CO re which we 
parametrise as: 

r core = r M 1/3 /(l + z) (5) 

Fig.|5]shows the cluster thermal SZ power spectrum for the 
gaussian (dotted line) and non-gaussian models (solid lines) , 
taking tp = 1 and r corc =0.08 h' 1 Mpc (for M = 10 15 ft _1 M o 
at z=0). The diamonds and vertical lines report the ACBAR 
and CBI observations, the star is the extrapolation of the 
primary CMB power spectrum, and the upper left curve 



is a fit to the WMAP measurements of the temperature 
power spectrum. In this plot, contrary to the previous para- 
graph, we did not require non-gaussian and gaussian models 
to have exactly the same abundance of massive clusters at 
2 = 0, but respect the normalisation given in Fig. The 
non-gaussian predictions clearly exceed the observations at 
£ ~ 2500, while the gaussian case underpredicts the power. 
To see how this depends on the free parameters of the mod- 
elling, Fig. |U] repeats the exercise for if> = and r core =0.15 
h~ x Mpc. The effect of changing these parameters is sig- 
nificant (mainly the increase in r corc ): predictions for the 
non-gaussian models agree with the observed power while 
the gaussian case has too little power with such a set of pa- 
rameters. At I — 2000, we note that the thermal SZ power 
spectrum shown in Fig.|U]for the g aussian case agrees within 
a factor 1.5 with that obtained bv lZhang et alJ JS)02) using 
hydrodynamical simulations. 

In fact, rather than the absolute level, it is the relative 
excess of the thermal SZ power spectrum in the non-gaussian 
models compared to the gaussian case which is interesting. 
Typically, at scales I J> 2000 non-gaussian models have more 
than an order of magnitude more power than the gaussian 
case. If as = 0.9 and if the excess power seen by CBI and 
ACBAR is mostly due to the cluster thermal SZ effect, this 
makes a strong case for non-gaussianity. Even if the excess 
factor we find for the power in the non-gaussian case were 
too high to fit observations given the local abundance of clus- 
ters, much room would be left for mildly non-gaussian mod- 
els with only reduced departure from gaussianity (the V and 
C 2 simulations used here have significant non-gaussianity). 
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Figure 8. WMAP temperature angular power spectrum (upper 
left curve) compared with the cluster thermal SZ power spectrum 
computed in non-gaussian simulations (solid lines for V and C 2 ) 
and gaussian (dotted line) models for clusters. The M — T relation 
assumes here ip — 1 and the /?-model r CO rc — 0.08 h,- 1 Mpc. The 
diamonds and vertical lines report the observations and the star is 
the extrapolation of the primary CMB power spectrum. Positive 
evolution in M — T and compact clusters increase the power at 
small scales. 
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Figure 9. Same as Fig. but for a M — T relation with -0 = 
and a /3-model with r C orc— 0.15 h- 1 Mpc. 

4 CONCLUSIONS 

While CMB-scale observations have started to constrain 
the possible amounts of non-gaussianity in the primor- 
dial cosmic density field, models with scale-dependent 
non-gaussianity could remain attractive for several as- 
pects of structure formation if they predict significant 
non-gaussianity on cluster scales. Here, we have addressed 
whether such non-gaussianity with reasonable scaling rela- 
tions can indeed explain the observables at high and low 
redshifts. 

We have confronted three recent observational hints of 
early cluster formation and reduced late-time X-ray clus- 
ter abundance evolution to the predictions of two phe- 



nomenological models with strong non-gaussianity on galaxy 
to supercluster scales and opposite amounts of primordial 
skewness. Our approach using collisionless simulations was 
mainly illustrative; for instance, we did not fit the parame- 
ters of the models we employed so that they precisely match 
the present-day galaxy power spectrum, but we have simply 
normalised the models so that they approximately reproduce 
the abundance of local clusters predicted in the concordance 
ACDM cosmology. 

Rather, we probed whether strongly non-gaussian mod- 
els such as a model with primordial voids with negatively 
skewed overdensity distribution, and a x-square model with 
one degree of freedom with positively skewed overdensity 
distribution would help resolve three possible difficulties 
faced by the gaussian ACDM paradigm. 

Firstly, we h ave found the protocluster observed by 
iMilev et alJ (I2004T) to be significantly easier to form in the 
two non-gaussian models we have simulated than in the 
ACDM case, where its estimated velocity dispersion seems 
difficult to reach. 

Secondly, we have compared the slow bright X-ray clus- 
ter abundance evolution observed in the 160 deg 2 ROSAT 
survey at z <J 0.6 to the predictions of the non-gaussian 
models using two models for the redshift evolution of the 
L — M relation. To reproduce the data, the simulated gaus- 
sian model required significant L — M evolution, which is 
observationally disfavoured. In the case of the two non- 
gaussian models, both evolving and non-evolving L — M 
relations provided predictions bracketing the data. 

Thirdly, we tackled the excess of power in CMB temper- 
ature maps seen by the interferometers at I ~ 2500. Using 
a simple /3-model for the intracluster gas distribution, we 
found that the thermal SZ contribution to the CMB tem- 
perature power spectrum is typically an order of magnitude 
higher in the non-gaussian models than in the gaussian case, 
irrespective of the two normalisations we have tried for the 
core radius of the intracluster gas and for the evolution of 
the T — M relation. If the observed excess of power is due 
to the thermal SZ effect, non-gaussianity may be required 
for ACDM to also match the local abundance of clusters (an 
alternative being an adjustment of the scaling relations). In 
fact, more "moderate" amounts of non-gaussianity than the 
ones considered here might suffice, and the room between 
the two strongly non-gaussian models we have simulated 
and a gaussian primordial density field should be explored 
in future work. 

To summarise, we have found that cluster-scale non- 
gaussianity can simultaneously explain an important con- 
straint at high redshift, and two observables at lower redshift 
when we use proper scaling relations. 

New constraints on possible non-gaussianity of the pri- 
mordial density field on relatively small scales will come 
from large weak lensing surveys. They are expected by 
lAmara fc Refreeierl i2005t) to yield constraints at a level sim- 
ilar to those set by current all-sky CMB maps on very large 
scales. In addition, both the inner density profiles of large 
DM haloes and the abundance of substructure are sensitive 
to the sign of the skewness of the prim ordial graviational 
potential field iAvila-Reese et all [2003'). Finally, new mea- 
surements of the 2 and n-point correlation function of local 
galaxy clusters of different masses will also set tight con- 
straints. If they yield negligible cluster scale non-gaussianity 
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then cluster formation and evolution might need substantial 
revisions. If not, simulations such as those presented here 
will need to be generalised. 
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